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Abstract 

The observation of charged lepton flavor non-conservation would be a clear signature 
of physics beyond the Standard Model. In particular, supersymmetric (SUSY) models 
introduce mixings in the sneutrino and the charged slepton sectors which could imply 
flavor-changing processes at rates accessible to upcoming experiments. In this paper we 
analyze the possibility to observe Z — > £j£j in the GigaZ option of TESLA at DESY. 
We show that although models with SUSY masses above the current limits could pre- 
dict a branching ratio BR(Z — ► fie) accessible to the experiment, they would imply an 
unobserved rate of [i — > ej and thus are excluded. In models with a small mixing angle 
between the first and the third (or the second and the third) slepton families GigaZ could 
observe Z — > tu (or Z — ► re) consistently with present bounds on lj — > £7-7. In contrast, 
if the mixing angles between the three slepton families are large the bounds from \i — ► e~f 
push these processes below the reach of GigaZ. We show that in this case the masses of 
the three slepton families must be strongly degenerated (with mass differences of order 
10 -3 ). We update the limits on the slepton mass insertions 5ll,rb,,lr and discuss the 
correlation between flavor changing and — 2 in SUSY models. 



PACS numbers: 12.60. Jv, 13.35.-r, 13.38.Dg 



1 Introduction 



Lepton flavor violation (LFV) has been searched in several experiments. The current status in 
\x and r decays is 
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In Z decays we have 



BR(Z->/ie) < 1.7 x 10~ 6 [6], 
BR(Z -> re) < 9.8 x 10~ 6 [6] , 

BR(Z -> r/i) < 1.2 x 10~ 5 [7] . (3) 

These observations are obviously in agreement with the Standard Model (SM), where lepton 
flavor number is (perturbatively) conserved. 

On the other hand, neutrino oscillations are a first evidence of LFV. Small neutrino masses 
and mixings of order one suggest the existence of a new scale around 10 12 GeV [8]. Massive 
neutrinos could be naturally accommodated within the SM (the so called z/SM). The contri- 
butions from the light neutrino sector to other LFV processes, however, would be very small: 
BK(£j — > ij'-f) < 10~ 48 and BR(Z — > liij) < 10~ 54 [9]. In consequence, any experimental 
signature of LFV in the charged sector would be a clear signature of nonstandard physics. 

In this paper we will study the implications of supersymmetry (SUSY) onZ^ P-iP-j- 1 The 
GigaZ option of the TESLA Linear Collider project [11] could reduce the LEP bounds down 
to [12] 
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with k = 0.2 — 1.0. We will here explore the possibility that SUSY provides a signal accessible 
to GigaZ in consistency with current bounds from BR(£j — * £77). Note that in SUSY models 
the branching ratio BR(^j — > 3£j) ~ a em BR(£j — > £77) will place weaker bounds on SUSY 
parameters (see Eqs. (1,2)) The conversion rate \i — > e on Ti gives also weaker bounds at 
current experiments, although this may change in the future (see [13] for a recent review). 

1 A recent work on the flavor-changing decays Z — > didj in 2HDMs and SUSY has been presented in [10]. 
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We will concentrate on the minimal SUSY extension of the SM (MSSM) with R-parity and 
general soft SUSY-breaking terms. Related works on LFV in Z decays in SUSY models study 
the MSSM [14] and a left-right SUSY model [15]. Several groups have analyzed other LFV 
processes in SUSY grand unified models with massive neutrinos (motivated by the atmospheric 
and solar neutrino anomalies [16]), or with R-parity violation [17]. There are also studies [18, 19] 
relating LFV Z decays with other processes. Direct signals of lepton flavor non-conservation 
in slepton production at the LHC [20] and at future e + e~ or fi + fi~ colliders [21] have been also 
explored. 

Other works on LFV Z decays in alternative models include the SM with massive Dirac or 
Majorana neutrinos [9, 22], left-right symmetric models [23], models with a heavy Z' boson 
[24], two Higgs doublet models (2HDMs) [25] or technicolor [26]. 



2 Calculation 



The most general vertex Viitj coupling a (lepton) fermion current to a vector boson can be 
parametrized in terms of four form factors: 



M = ie^u £l (p 2 ) [l^{Fv - F A >y 5 ) + (iF M + F E ^)a liV q v \ u tj {j)\) 



(5) 



where Ey is the polarization vector (ey • q — 0) and q = p 2 — p± is the momentum transfer. 
For an on-shell (massless) photon Fa = 0, and, in addition, if mj ^ mj then Fy = 0. This 
implies that the flavor-changing process £j — > £77 is determined by (chirality flipping) dipole 
transitions only. In contrast, all form factors contribute to the decay of a Z boson: 



BR(Z->^j) = BR(Z -> lilj) + BR(Z -> ijij) 
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with a w = g 2 /(4:7i), f LtR = f v ±fA, fv,A = -I6n 2 g 3 F VtA and f M ,E = -16M w n 2 g 
calculate (see Appendices A and B for details) these branching ratios in the MSSM 
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We 



Let us consider the case with two lepton families. Since SUSY is broken, fermion and 
scalar mass matrices will be diagonalized by different rotations in flavor space. After the 
diagonalization of the fermion sector we are left with a 2 x 2 scalar matrix with 3 arbitrary 
parameters. We will assume that the rotation that diagonalizes the scalar matrix is maximal, 
6 = 7r/4 (i.e. we assume no alignment between fermion and scalar fields; the amplitudes that 
we will calculate are proportional to sin 29). Our choice corresponds to a mass matrix with 
identical diagonal terms. Taking 



2 ~ 2 
m = m 



5 



5 



the two mass eigenvalues are 



m 2 



1,2 



= m 2 (v / TT^T5) • 
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In this parametrization m 2 = mim 2 characterizes the SUSY-breaking scale and 5 = (m| — 
ml) / (2m 2 ) the mass splitting between the two families. 5 is also responsible for any flavor- 
changing process: 5 = corresponds to the flavor-conserving case, 5 <C 1 can be treated as a 
non-diagonal mass insertion, and 5 — > oo gives m\ — > oo (a decoupled second family). The last 
case implies a maximum flavor-changing rate [14, 15]. 

To analyze the general case with three lepton families we will consider two scenarios. First, 
we will follow the usual approach [27] where the influence of a nondiagonal term 5 IJ is calculated 
putting the rest to zero. This implies that the slepton family £k^i,j does not mix with £j and 
£j (Qik = Qjk = 0), which reduces the problem to the two family case discussed above. This 
approximation is only justified if the off-diagonal terms satisfy 5 IK 5 KJ < 5 IJ or, in terms of 
mixing angles and mass differences, if 

f mj-m 2 K \ ( m\-m 2 \ m 2 - m 2 

~ o I - 2 S111 U IK r~^T S111 U JK ~ - 2 , - 2 S111 

ymj+m^ J \m K + mj J mj + mj 

We will then discuss a second scenario with maximal mixing between the three slepton families: 
#12 = #23 = #13 = 7r /4. Large mixings are suggested by the observation of solar and atmospheric 
neutrino oscilations (note, however, that the non observation of v e — v T oscilations in CHOOZ 
[28] could suggest 9 n < 0.1). We will use 5 IJ to parametrize the mass difference between i T 
and if. 5 IJ = (m 2 — m 2 )/(2m\), where m\ is the mass of the lightest slepton family. 

The relevant parameters for the calculation will then be the masses and mixings of charginos 
and neutralinos, the masses of the six ('left' and 'right' handed) charged sleptons, and the 
masses of the three sneutrinos. When we evaluate the contribution of each 5 IJ to Z — > ijlj 
and £j — > ij'-f setting all the other to zero we will have three independent parameters 5 1 [ I L J 
in the sneutrino sector and nine 5^ , 5^ and 5 L ^ in the charged slepton sector. In the 
case of maximal mixing between the three slepton families there will be two independent mass 
differences <5^ J and four 5 1 / 1 / 7 , 5 e R ^ (note that in this case 5 23 = 5 13 — 5 12 ). 

In our analysis we will not assume any (grand unification) relation between slepton masses. 
For each non-zero choice of 5 IJ it is straightforward to obtain and diagonalize the mass matrix 
that corresponds to a maximal rotation angle (see Appendix B). Our results should coincide 
with the ones obtained in the limit of small mass difference using the mass insertion method, 
but they are also valid for any large value of 5 IJ . 

The process Z — > £j£j goes through the diagrams in Fig. 1. Box diagrams mediating 
e + e~ — * ill j introduce a small correction of order V z/Mz ~ a. Analogous diagrams describe 
£j — > ij'y. The inclusion of the contributions of the third type is essential to cancel ultraviolet 
divergences (they are related to counterterms by Ward identities). Diagrams with neutralinos 
in (A) or sneutrinos in (B) do not couple to the photon. The diagrams of type (C) do not give 
dipole contributions. 

Due to the weaker experimental bounds (in Table 1) on sneutrino masses, the dominant 
contributions to Z — > tjlj will come from the diagrams mediated by chargino-sneutrino (see 
Fig. 2). Note that sneutrino masses can be substantially lighter than charged slepton masses 
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Figure 2: BR(Z — + fie) (lower curves) and BR(yU — > ej) (upper curves) as a function of the 
lightest scalar mass rh\ for tan/5 = 2 and the different S 12 . Solid lines correspond to M2 = 150 
GeV and /i = —500 GeV and dashed lines to M 2 = /i = 150 GeV. 
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Table 1: Approximate lower bounds on SUSY mass parameters based on [29]. Note that, for 
negligible scalar trilinears, to? = to- + M§ c^y cos 2/3, and the bounds on mp and to^ are 
correlated. For instance: rxiy > 65 (40) GeV for tan/3 = 2 (50). 



lightest slepton (mi) 


m 9 > 45 GeV 
m/ > 90 GeV 


lightest chargino 


> 75 GeV, if rrii, > m + 

> 45 GeV, otherwise 


lightest neutralino 


m^o > 35 GeV 

Xi — 



for large tan (3 and light SUSY-breaking masses, 

to- « 
to| « 

which tends to increase the maximum relative contribution of chargino-sneutrino diagrams. 

We would like to emphasize that our results will depend on contributions with opposite 
signs that often cancel when varying a parameter. For example, one would expect that the 
process Z — > Ijlj is optimized for light slepton masses. However, we observe frequently the 
opposite effect. Its branching ratio can increase by raising the mass of the sleptons up to values 
of 500 GeV, and only at masses above 1 — 2 TeV the asymptotic regime is reached (see Fig. 2). 
These cancellations give a one or two orders of magnitude uncertainty to any naive estimate, 
and underline the need for a complete calculation like the one presented here. 

We give in Fig. 3 the dominant diagrams in terms of gauginos, current eigenstates and mass 
insertions, specifying the chirality of the external fermion. All the diagrams contributing to 
£j — > £/7 except for the last one grow with tan (5. 

3 Results 

3.1 Z -> tit j at TESLA GigaZ 

Let us consider the process Z — > iilj uncorrelated from other LFV processes. For SUSY 
masses above the current limits it is possible to have Z — > fie; re; t/j, at the reach of GigaZ. 
The maximun rate is obtained when the second slepton £j is very heavy (i.e. 5 IJ — > oo). The 
largest contribution comes from virtual sneutrino-chargino diagrams (all other contributions 
are at least one order of magnitude smaller). It gives BR(Z — > ijtj) from 2.5 x 10 -8 for 



m 2 L + -M 2 Z cos 2(3 , 



»>l + \ -- + s 2 w ) M| cos 2/3, 



(10) 
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Figure 3: Dominant diagrams contributing to (a) Z — > fie and (b) fx —>■ erf, in terms of 
gauginos, higgsinos and current eigenstates, showing the approximate linear dependence on the 
flavor-changing mass insertions 5 12 (crosses), the fermion mass insertions (big dots) and tan/3. 

tan/3 = 2 to 7.5 x 10~ 8 for tan/3 = 50, practically independent of the lepton masses. The 
variation is due to the mild dependence of chargino and sneutrino masses on tan (3. These 
branching ratios are above the values given in Eq. (4). We find that a branching ratio larger 
than 2 x 10~ 9 (2 x 10~ 8 ) can be obtained with sneutrino masses of up to 305 GeV (85 GeV) 
and chargino masses of up to 270 GeV (105 GeV). 

Most of these values of BR(Z — > £i£j), however, are correlated with an experimentally 
excluded rate of £j — > £7-7. We give below the results in the two scenarios (independent off- 
diagonal terms and maximal mixing of the three flavors) described in the previous section. 

(%) We separate the contribution of each 5 IJ setting all the other to zero. For the first 
two families, after scanning for all the parameters in the model we find that BR(/i — > e^) < 
1.2 x 10 -11 implies BR(Z — > fie) < 1.5 x 10~ 10 , which is below the reach of GigaZ. 

A more promising result is obtained for the processes involving the r lepton. It turns out 
(see also next section) that the bounds from r — > £7-7 can be avoided while still keeping a rate 
of Z — ► t£j at the reach of the best GigaZ projection (see Fig. 4). In particular, for large 
5ll 3 (or 5 V L 2 £") and a light sneutrino (of around 70 GeV) we get BR(Z — > re) w 1.6 x 10~ 8 
for BR(r — > ej) ~ 3.5 x 10~ 8 , which is two orders of magnitude below current limits (with 
similar results for BR(Z — > Tfi) and BR(r — > ffy)). This result is due to the sneutrino-chargino 
diagram. The contributions due to charged slepton mixing are essentially different in the sense 
that they saturate the experimental bound to r — > £j^ giving a small effect (at most, one order 
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Figure 4: BR(Z — > t/i) and BR{r — > 717) as a function of the lightest sneutrino mass (mi) 
with the other one decoupled (5£f 3 — > 00), in several SUSY scenarios at the reach of GigaZ. 

of magnitude below the reach of GigaZ) in Z — > t£j. We obtain events at the reach of GigaZ 
with lightest sneutrino masses from 55 to 215 GeV, lightest chargino from 75 to 100 GeV, and 
tan (3 up to 7. 

(ii) In the case with maximal mixing between the three slepton flavors it is not consistent 
to take 5 IJ 7^ and S IK = 5 KJ = 0. In terms of slepton mass differences, only two of the 
three mass differences are independent (5 23 = 5 13 — 5 12 ). In terms of off-diagonal terms 5 IJ 
in the mass matrix, for maximal mixing only one of them can be put to zero. Note that in 
this case the non-observation of // — ■> e'-f will constraint all the 5 IJ parameters, not only 5 12 : a 
non-diagonal S 12 mass insertion would be generated through a 5 13 followed by a 5 32 . In fact, 
we find that the constraints from r — > e'j; /J7 are always weaker than the one from fi e^y. A 
branching ratio BR(/i — > cy) < 1.2 x 10 -11 implies BR(Z — > fie; re; t/j) < 10 -9 , and the three 
lepton flavor violating decays of the Z boson would be out of the reach of Giga Z. 



3.2 Bounds on 5 IJ from £j — > £77 

The bounds on the <5 /J parameters establish how severe is the flavor problem in the lepton 
sector of the MSSM. We will update them here, including the sneutrino-chargino contributions 
neglected in previous works [27] and the general slepton-neutralino contributions (photino 
diagrams are typically subdominant as pointed out by Ref. [30]). In addition, we also consider 
the case of maximum mixing between the three slepton families. 

The limits come exclusively from the process ij — > £77. To estimate the MSSM prediction 
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Table 2: Bounds on the 5 12 s from BR(/i — > e7) < 1.2 x 10~ n in different SUSY scenarios 
assuming no mixing with the third family. The bounds on 5 IJ in the case of three family 
mixing can be read from these ones (see text). 



Svl2 
LL 


mi 


M 2 


li = -500 


H = -150 


fx = 150 


LI = 500 


tan = 2 


100 


150 


14 x 10~ 3 


1.0 x 10~ 3 


0.3 x 10~ 3 


1.0 x 10~ 3 


500 


33 x 10" 3 


3.0 x 10~ 3 


1.7 x 10~ 3 


13 x 10" 3 


500 


150 


3.7 x 10~ 3 


1.1 x 10~ 3 


1.2 x 10~ 3 


1.7 x 10~ 3 


500 


7.3 x 10~ 3 


2.6 x 10~ 3 


2.3 x 10~ 3 


4.1 x 10~ 3 


tan p = 50 


100 


150 


9.3 x 10~ 5 


2.1 x 10~ 5 


2.0 x 10~ 5 


8.6 x 10~ 5 


500 


80 x 10~ b 


9.0 x 10~ 5 


8.8 x 10- 5 


77 x 10" b 


500 


150 


9.7 x 10~ b 


4.5 x 10- & 


4.5 x 10~ b 


9.4 x 10~ b 


500 


22 x 10" b 


9.6 x 10" b 


9.5 x 10~ b 


21 x 10" b 


LL 


rh\ 


M 2 


li = -500 


fx =- -150 


li = 150 


H = 500 


tan = 2 


100 


150 


7.5 x 10~ 3 


2.6 x 10~ 3 


1.7 x 10~ 3 


5.0 x 10~ 3 


500 


84 x 10~ 3 


11 x 10" 3 


7.0 x 10~ 3 


41 x 10" 3 


500 


150 


14 x 10~ 3 


8.5 x 10" 3 


0.12 


32 x 10~ 3 


500 


ij — 

24 x 10~ 3 


T 

20 x 10^ 3 


n 

32 x 10~ 3 


n 

26 x 10~ 3 


tan p = 50 


100 


150 


2.4 x 10~ 4 


0.8 x 10" 4 


0.8 x 10~ 4 


2.4 x 10~ 4 


500 


2 

24 x 10" 4 


1 — 

3.5 x 10" 4 


T 

3.4 x 10" 4 


T 

23 x 10" 4 


500 


150 


— — 1 — 

7.4 x 10~ 4 


: 1 — 

6.4 x 10" 4 


T 

6.9 x 10~ 4 


— — 1 — 

7.6 x 10" 4 


500 


9.7 x 10 -4 


9.6 x 10~ 4 


9.8 x 10~ 4 


9.7 x 10~ 4 


c/ 1 1 9 


rhi 


M 2 


ix = -500 


li = -150 


li = 150 


H = 500 


tan /3 = 2 


100 


150 


4.2 x 10~ 3 


1.5 x 10~ 3 


1.8 x 10~ 3 


3.7 x 10~ 3 


500 


11 x 10" 3 


3.2 x 10~ 3 


2.5 x 10~ 3 


8.3 x 10~ 3 


500 


150 


22 x 10" 3 


10 x 10" 3 


22 x 10" 3 


0.33 


500 


19 x 10" 3 


12 x 10" 3 


0.33 


35 x 10" 3 


tan /3 = 50 


100 


150 


1.6 x 10~ 4 


0.6 x 10" 4 


0.6 x 10~ 4 


1.5 x 10~ 4 


500 


3.8 x 10" 4 


1.1 x 10" 4 


1.1 x 10" 4 


3.8 x 10~ 4 


500 


150 


1 — 

16 x 10~ 4 


1 1 

13 x 10" 4 


7\ — 

15 x 10" 4 


— 1 — 

17 x 10~ 4 


500 


9.4 x 10~ 4 


8.1 x 10~ 4 


8.7 x 10~ 4 


9.6 x 10~ 4 


LR 


rhi 


M 2 


li = -500 


ix = -150 


li = 150 


li = 500 


tan /3 = 2 


100 


150 


1.6 x 10~ 6 


1.5 x 10~ 6 


1.6 x 10- 6 


1.7 x 10~ 6 


500 


4.5 x 10-° 


4.4 x 10~ d 


4.7 x 10-° 


4.6 x 10~ B 


500 


150 


1.3 x 10-° 


1.2 x 10~ d 


1.2 x 10-° 


1.2 x 10~ B 




500 


7.6 x 10-° 


7.5 x 10~ d 


7.6 x 10-° 


7.7 x 10~ B 


tan /3 = 50 


100 


150 


1.6 x 10- 6 


1.5 x 10~ 6 


1.6 x 10~ 6 


1.6 x 10" 6 


500 


4.5 x 10~ s 


4.5 x 10~ d 


4.6 x 10-° 


4.5 x 10~ B 


500 


150 


1.3 x 10~ e 


1.2 x 10- y 


1.2 x 10~ B 


1.3 x 10~ e 


500 


7.7 x 10~ e 


7.6 x 10~ d 


7.6 x 10~ e 


7.7 x 10- e 



8 



we combine low and high values of the relevant parameters: tan (3 — 2; 50, mi = 100; 500 GeV, 
and the gaugino and higgsino mass parameters M 2 = 150; 500 GeV and \x = ±150; ±500 GeV. 

(i) The results for the case with a decoupled family are summarized in Table 2. We include 
the bounds from /i — > to S^i, 5 RR and S^r. A 5 ~ 10~ 3 implies a l%o degeneracy 
between the two slepton masses. We observe that the degeneracy between the selectron and 
the smuon is required even for large SUSY masses, and it must be stronger if tan/3 is large, 
as expected from the diagrams in Fig. 3. The small values of Sj} R , around 10 -6 , imply just 
that the scalar trilinears, usually assumed proportional to the Yukawa couplings, are small. 
Particularly weak bounds on the <5's (bold faced in Table 2) are obtained when approaching 
a dip of the curves in Fig. 2. This occurs for certain values of the SUSY parameters due to 
cancellations of the contributions of the various particles running in the loops. 

The experimental bounds on the mass differences involving the third family are much weaker. 
They come from r — > £77 (and not from /i — > e^f), since in this case we are asuming that the 
mixing with £j (J 7^ /, 3) is negligible. For small tan (3 we find no bounds on any 5 13 (except for 
8 LR ). For large tan/3 the bounds are (depending on the values of the SUSY-breaking masses) 
5££ 3 = 0.03 to 1.3; <5£l = °- 14 to °°; and 5 rr = °- n to 00 ■ For the LR mass insertions, 
we find 8 LR = 0.05 to 00, independent of tan/3. This results improve the bounds obtained in 
Ref. [27], in particular, the ones involving the second and third lepton families. 

(ii) As explained before, if the three slepton families are maximally mixed (as suggested 
by the experiments on neutrino oscilations) then the strongest bounds on all the slepton mass 
differences come from /i — > exclusively. As we will see, however, the bounds can be read 
directly from Table 2, since the small values obtained for 5 12 admit an analysis based on mass 
insertions. Let us first suppose that the first and second sneutrino families are degenerated 
(^lz, 2 = 0)- Then a mass difference 5^ = 51 23 = 5 generates an off-diagonal mass term 
(see Appendix B3) S^} 2 = 5/4. In consecuence, the degree of mass degeneracy with the third 
family (5) imposed by /i — > is just four times the value given in Table 2. We have checked 
numerically that this estimate is quite accurate. Analogously, the bounds on = 8 V L ]^ when 
^ll 3 = an d on b v j} 2 = (^f 3 when 5^ = are respectively two and four times the values in 
Table 2. In the same way we can read there the bounds on 5 e L J £ J and 5 R T ^, which establish the 
degree of mass degeneracy between the three families of charged sleptons. 

3.3 Lepton flavor violation and g M — 2 

Finally we would like to comment on the relation between /j, — > e^y and the muon anomalous 
magnetic dipole moment. See Ref. [31] for more exhaustive analyses of the constraints on lepton 
flavor violation in the MSSM from the muon anomalous magnetic moment measurement. A 
<7 M — 2 correction would be generated by the diagrams in Fig. 3b if no mass insertions 5^£ J , 
^lli Srr are included and 5 e LR J is replaced by S e L 2 R . In this sense, — 2 is a normalization of 
the branching ratio BR(£j — > £77) for processes changing the muon flavor. 

We plot in Fig. 5 the value of a M = (g^—2)/2 for the SUSY parameters in the region accessible 
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Figure 5: Total SUSY contribution to the muon anomalous magnetic dipole moment as a 
function of the sneutrino mass with m L = m R , in the same SUSY scenarios as in Fig. 4. 

to GigaZ not excluded by r — > /ry, taking for simplicity equal soft-breaking terms mj, = rriR 
(they would not very different, for example, assuming left-right unification at the GUT scale). 
We obtain, in agreement with [32], positive or negative contributions correlated with the sign of 
the Higgsino mass parameter \x and similar in size to the weak corrections. The recently revised 
SM prediction [33], a^ M = 11 659 179.2 (9.4) x 1CT 10 , compared to the world average after the 
last data from the Brookhaven E821 experiment [34], a c * p = 11 659 202.3 (15.1) x 10~ 10 , exhibits 
a 1.4cr discrepancy: 5a^ = a° xp - a^ M = (23.1 ± 16.9) x 10~ 10 . This indicates that the muon 
dipole moment may still need non-standard contributions of positive sign. In any case, the 
MSSM contribution <5a^ USY is bounded at two standard deviations by the dotted lines in Fig. 5. 
Only the regions with heavier masses in the scenarios of Fig. 4 are favored. 

Fig. 6 shows the muon anomalous magnetic moment for the different sets of SUSY param- 
eters employed to explore the muon LFV decay in Table 2. Low values of tan/5 and positive 
values of ji are preferred by — 2, which implies less stringent bounds on the 5 insertions 
parametrizing the flavor-changing lepton decays. 



4 Conclusions 

SUSY models introduce LFV corrections which are proportional to slepton mass squared dif- 
ferences. We have shown that the non-observation of fi — > e7 implies around a l%c degeneracy 
between the masses of the sleptons in the first two families. Once this degeneracy is imposed, 
the rate of Z — > fie is always below the limits to be explored at GigaZ. Moreover, if the mixing 
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Figure 6: Total SUSY contribution to the muon anomalous magnetic dipole moment as a 
function of the sneutrino mass, with uil = m,R. 

between the three slepton families is large then also the third family must be (at the l%o) 
degenerated, and the processes Z — > re; r/i will not be observed at the GigaZ. The degeneracy 
between the lightest slepton families could be justified by the weakness of its Yukawa couplings, 
but for the third family it should put constraints on definite SUSY models. 

In contrast, if the mixing between the third and the first slepton families is small, then the 
third family could be much heavier than the other two and there would be no flavor problem in 
the slepton sector (the bounds would come only from r — > /ry, not from /i — * cy). In this case, 
if the GigaZ option of TESLA reaches its best projected sensitivity it could observe Z — > t/i 
coming from the virtual exchange of wino-sneutrino. 
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Note added 

After the completion of this paper, new data from BNL g — 2 experiment [41] appeared con- 
firming the previous measurement with twice the precision. The discrepancy with the SM 
prediction is now more significant, up to 2.6a. There are also new experimental bounds on 
t — > /ry [42]. The new data do not introduce any qualitative changes in our results. 
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A Generic expressions at one loop for Zijij 

A.l Feynman rules in terms of generic vertex couplings 

Let / be a fermion, a scalar field and Pr,l — (1 ± 7s)/2- The Feynman rules for the three 
vertex topologies needed are: 



• Vertex Zf A f B 

• Vertex Z4> x (p 2 )<Py{pi) 

• Vertex 4>xfAfi 



^gi^g L ABPL + g%P R ) , 

igGxviPi +P2Y , 
W^axPl + cf AX P R ) . 



A. 2 Invariant amplitude 



The most general invariant amplitude for on-shell external legs is 



Mi 



w 



Let us introduce the squared mass ratios \ n = m^/M^ and the dimensionless two- and three- 
point one-loop integrals 



Bi(A ,Ai) = B 1 (0;m 2 ,m 2 1 ) , 
C..(A ,Ai,A 2 ) = M$y C..(0, M§, 0; ml, m\, ml) , 



(12) 
(13) 



from the usual tensor integrals [35, 36], 



B»(3?-mlm\) = p»B x , (14) 
C»(p 2 1 ,q 2 ,p 2 2 ;m 2 ,m 2 1 ,m 2 2 ) = tfC n +l%C 12 , (15) 
C^(p?,<^pl;mg,m?,™l) = p ^C 21 +p^C 22 + (rftf + p^)C 23 + g^M^C u . (16) 

Note that C , C23, C 24 , Cn + C i2 and C 2 i +C 22 are symmetric under the replacements Ai <-> A 2 , 
while Cn — C i2 and C 2 i — C 22 are antisymmetric. The form factors for each type of diagrams 
(Fig. 1) are: 

• Diagram of type A: 



fl XS = E [vab^ax^bx V / A^C (Ax ! A a ,Ab) 

XAB ^ 



Az C 2 3(Ax, Aa, A#) — 2 C 2 4(Ax, A^, \b) 



fXXS = f XXS {L ^ Rh 



(17) 
(18) 
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/m* — E { 2 [(9AB C IAX C JBX + 9AB c fAX c JBx) + (I <-> J)*] 
XAB ^~ 



x [C n (A x , A A , A b ) + C 21 (Ax, A A , A B ) + C 23 (Ax, A A , Ab)] 



+^ [(9abCiaxC% x + 9abC?axC% x ) + <- •/)*] Ci 2 (Ax, A A , A B ) 



/s* — E { 2 ^9ab c iax c jbx 9ab c ?ax cI jbx) (I J)* 

XAB ^~ 



x [Cn(Ax, A A , Ab) + C 2 i(Ax, A A , Ab) + C 23 (Ax, Xa, Ab) 
A7 



K^abC/axCjbx " 9AB<$Axtj B x) -(/<-> J)*] C 12 (A X , A A , A B ) 

Diagram of type B: 

fL X = ~ 2 G X yCiax c jay ^24(Aa, Ax, Ay) , 
A xy 

f s M x = E { - ^ Gxy ^ cL ^ + c *a* c *ay) + (i <- jy\ 

AXY ^ 

x [Cn(A A , A x , Ay) + C 12 (A A , Ax, Ay) + C 23 (A A , Ax, Ay)] 
£± [Gxyc^ x cV + (/ - J)*] 

x [C (A A , Ax, Ay) + C n (A A , Ax, Ay) + C 12 (A A , A x , Ay)] | , 

if s E x = { ~ ^ ^y^xc L jAY - cf: x c^ AY ) - (I ^ J)*] 

A YV V 

[C n (A A , Ax, Ay) + C 12 (A A , Ax, Ay) + C 23 (A A , Ax, Ay)] 

^ [G X ycY AX C* AY - (/ - J)*] 

[C (A A , Ax, Ay) + Cn(A A , Ax, Ay) + Ci 2 (A A , Ax, Ay 



AXY 



2 

< 



Diagram of type C: 



ff = - Q ^Y. C IAXC L JAXM^XX) , 

s 2 

/if = — / ; C iax c jax Bi(A A ,Ax) , 
Cw 7t 

f X M = 0, 

ff = o. 
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The tensor integrals are numerically evaluated with the computer program LoopTools [37], 
based on FF [38]. 

Non-trivial checks of our expressions are the finiteness of the amplitude and the test of 
the decoupling of heavy particles running in the loops, that must take place both in the SM 
and the MSSM [39]. These conditions are fulfilled only when summing over the different type 
of diagrams involved thanks to the relations existing among vertex couplings. Note that the 
ultraviolet-divergent tensor integrals are the same that diverge with a large mass M, 

C 24 -> ~ - ^logM, Bi^- + logM, e = D-A. (29) 

All the other tensor integrals are finite and vanish for large masses. 



B Masses, mixings and vertex couplings in the MSSM 

Notation: the indices I or J refer to the flavor of the external fermion; the indices A or B refer 
to a chargino/neutralino mass eigenstate (xa=i2 an< ^ Xa=i 2 3 4)> the indices X or Y refer to a 
charged slepton/sneutrino mass eigenstate (£x=i,...,6 an d ^=1,2,3) ■ 



B.l Charged sleptons 



Let £l z and £ Rl be the superpartners of the charged leptons £l 2 and £ Rl , respectively. The 6x6 
mass matrix of three generations of (charged) sleptons can be written as 



Mf 



m 



LL 



m 



2 T 
LR 



m L m RR 



(30) 



where m| L and m% R are 3x3 hermitian matrices and m LR is a 3 x 3 matrix, given by 



(m£ L ) J<7 = (ml)! j + ^71% + \ -- + s 2 w j M 2 Z cos 2(3 
( m D/J = i m R)u + \ m \ - M z cos2l3s 2 w ]5 u , 

/ 2 \ 1 * \ v C0S P a r 

( m Li?)/J = (A*)jj— -j= mtj/itaap d u . 

The mass matrix M| can be diagonalized by a 6 x 6 unitary matrix S^, 

S £ M|S ft = diag(m| x ), X = l,...,6. 
The mass eigenstates are then given by 

^ = Si/ L/ + , X = 1, . . . , 6, / = 1, 2, 3 . 



(31) 
(32) 
(33) 



(34) 



(35) 
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B.2 Sneutrinos 



There are only 'left-handed' sneutrinos in the MSSM. Let u Ll be the superpartner of the 
left-handed neutrino uj. Then the 3x3 sneutrino mass matrix contains the same soft SUSY- 
breaking mass term as the 'left-handed' sleptons and a different D term: 

(MDj j = (ml) j j + \m\ cos 2(3 5 U , (36) 

and it is diagonalized by a 3 x 3 unitary matrix S^, 

M? = diag(m?J, X = 1,2,3, (37) 
so that the sneutrino mass eigenstates are 

vx = S'xjVLj, X = 1, 2, 3, I = 1, 2, 3 . (38) 



B.3 Slepton matrices in terms of 5 mass insertions 



(i) Assuming that only two generations (/ and J) of charged sleptons mix and they do it 
maximally (9 = 7r/4), only the following 4x4 symmetric mass matrix, with entries I,J,I + 
3, J + 3, is relevant: 



/ 1 



z,2 



Mj = m 



\ 



^ 1 " 



(39) 



The insertions 5 e L j^ and are flavor conserving. We assume that, alternatively, only one of 
these <5's is different from zero. Then, the relevant non-diagonal 2x2 submatrix: 



2 ~ 2 

m = m 



1 5 
5 1 



(40) 



is trivially diagonalized by the following submatrix of S: 



U 



-1 



(41) 



y/2 \ 1 1 
yielding the eigenvalues: 

m\ 2 = m 2 (Vl + S 2 =F 5) , (42) 
where 5 = (m| — mf) /(2m 2 ) is the mass splitting between both generations of sleptons. 
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The relevant 2x2 submatrix for the sneutrinos in terms of the mass insertion 8^[ J is 
constructed in a similar way. 

(ii) For the case when the three generations of sleptons mix we employ the standard 
parametrization for the relevant 3x3 submatrix: one CP phase (that we set to zero) and three 
mixing angles (9 12 , #13, #23 where 6u represents the mixing between familiess / and J when the 
mixing to the remaining one is zero. We take again maximal mixing, (9 12 = #13 = #23 = 7r /4. 
Then, 

( V2 -y/2 -2 \ 
y/2-1 V2 + 1 -y/2 
V2 + 1 V2-1 V2 



U 



2V2 



\ 



I 



is the unitary matrix that diagonalizes the symmetric mass matrix: 



m 



gl2 ~ 5 22 . 
5 13 P 3 ~5 33 



with 



(43) 



(44) 



6 11 


= l + V 2 + <5 13 )- 
4 


^(5 12 -5 13 ) 
4 


(45) 


5 22 


= l + V 2 + ( p) + 
4 


4 


(46) 




= l + i(5 12 + 5 13 ) 




(47) 


~5 12 






(48) 


~5 13 


= -\[(2-V2)5' 2 - 


(2 + V2)5 13 } 


(49) 


~ 6 23 


= -\[(2 + V2)5 12 - 


(2 - V2)5 13 ] 


(50) 


m 2 2 


= m 2 (l + 25 12 ), m\ 


= m 2 (l + 25 13 ). 


(51) 



yielding the eigenvalues: 



Now the mass splittings S IJ = (m 2 — m 2 )/(2m 2 ) are not the same as the off-diagonal mass 
insertions. 



B.4 Charginos 

The chargino mass matrix, in the (charged wino, charged Higgsino) basis, is 

X= I Ma V2M w smf3 
/% ~ " 1 V2M w cos(3 n 
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It can be diagonalized by two unitary matrices U and V, 



(53) 



where 



m 2 -± 

Xl,2 



Mj + /i 2 + 2M W 



=F yj (M| - /i 2 ) 2 + 4M W cos 2 2/5 + 4M 2 ,(M 2 + /i 2 + 2M 2 /i sin 2/3) 
In order to get positive-mass eigenstates, one introduces two orthogonal matrices O-t, 

U = o_ , 

0+ , detX > , 
cr 3 + , detX < , 

where a 3 is the usual Pauli matrix. 



V 



(54) 

(55) 
(56) 



B.5 Neutralinos 



(57) 



The neutralino mass matrix, in the basis of the U(l) and SU(2) neutral gauginos and the two 
neutral Higgsinos (B, W 3: Hf, H$), is the symmetric matrix: 

( M 1 ■ ■ ■ \ 

M 2 • • 

—MzSw cos (3 MzCw cos f3 

y M z s w sin (3 — M z c w sin P —fx J 

To simplify, we employ the unification constraint M\ = |tan 2 w M 2 . 

The matrix above can be numerically diagonalized by the unitary matrix N, 

N*YN _1 = diag(m^o, m^o, m^o, m^o) . (5£ 



B.6 Vertex couplings 

• Vertex ig^(g L AB P L + 9ab p r) [ note that 9ba = 9ab> 9ba = 9ab\- 

2 
1 



z XaXb 


9ab = 


1 n' L ■ 

u ABi 

cw 







9ab = 


—0 ,R ■ 

„ w ABi 
C W 





ZXaXb 


9ab = 


—o" L ■ 

U AB> 

C W 







9ab = 


1 o" R - 

u ABi 






IL 

AB 



- S w I \Ja2^J B 2 + C^UaiU^! 



iR 

AB 



llL 



nR 



— S 



W 



V^ 2 V B2 + c^V* 41 V 



W v Al v Bl 



1 

2 



'B3J ) 



y/L* 
AB 



(59) 
(60) 
(61) 
(62) 
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Vertex igGxy(pi + P2Y [note that Gyx = G* XY ] 

Zv x v Y : G X y = -77— $XY , 

2,Cw 



3 



1 _ « 2 \ at _ c 2 ce at 

2 b W ) °XK°YK b W J X,K+3 Y,K+3 



ZPJy ■ G XY = — L 
Vertex ig{c L lAX P L + cf AX P R ): 

^xXa^I : C IAX = ~ ^Al^XI > 

M% ■■ c 7 L JS = -^(tan^N^ + N^ 2 ) S XI - ™'* N^ 3 S^ +3 : 
V 2 v 2M W cos p 

y2M w cos p 



C The LFV decay £j £ n and p - 2 



The general amplitude £j — > £77 at one loop reads 

= -ie^e^^)— [{if M + fl^a^q"} u tj {px) 



47T 7 77l£ 



J 



In the literature one finds often the notation: 



,2 



For equal leptons, the anomalous magnetic dipole moment of I is 



The width of £j -> £77 is 



G^pffVg nctw 
Since the width — > tiVjvi) = ± and Gp = — — — — , one has 



192tt 3 " v^M, 2 



w 



BR(£;^77) _ 12a _ 48^ L 2 fl 2 

"U/mI + I /si J — ~~ — UA2 1 +nl ii 



BR(£j — > tjvjvi) 71 ml KlJMl lJhl ' G\ 
where BR(£j — > tiVjVj) = 1/0.17/0.17 for Ijlj = /le/refT/i, respectively. 
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The SUSY contributions to the form factors are the following. 
Diagram of type A: Chargino-Chargino-Sneutrino [xax = w^±/m~ ] 



f. 



mi, 



■ fl 



m £j 



vx 



A=l X=l 
2 3 

EE 



A=l X=l 



"Kj L[C]* L[C] „ ( v m xt L[C]* R[C] p / \ r p 
m u x Ul vx 



Diagram of type B: Slepton-Slepton-Neutralino [x° AX = m| l m \ x \- 



f7 
JM 



m tj 
where 



EE 

A=l X=l 
4 6 

EE 

A=l X=l 

F 2 (x) 
F 3 (x) 
F 4 (x) 



™ij L[N]* L[N] 



c^xF^Ax) + ?cSc;fjF 4 (4) + L ~ R 



rn 



^IAX ^JAX 



tx 



9 ^IAX ^.JAX J 
m ix 



m tj L[N]* L[N] p , o \ I "M L[N] * R[N] „ , n x _ r d 



771; 



2 + 3x — 6a; 2 + a; 3 + 6x In a; 

6(1 -x) A ' 
-3 + Ax - x 2 - 2 In x 
2(1 -x) 3 ' 
1 — 6x + 3a; 2 + 2a; 3 — 6x 2 In x 



6(1 -xf 
1 — x 2 + 2x In a; 



-xFi(l/x) , 



2(1 -a;) 3 

These functions are combinations of 3-point tensor integrals, in agreement with [40]: 

F 1 (x A x)/ml x = 2 [Cii + C 2 i + C 23 ] (0,0,0; m„ x ,m x ±, m^±) , 
,2 



F 2 (x AX )/m Px = 2 Cn(0, 0,0; m Px , m-±,m-±) , 
F 3 (a; AX )/m| x = -2 [C n + C 21 + C 23 ](0, 0, 0; m^m ix ,m ix ) , 
^(^x)/^ = [Co + C u + Ci 2 ](0, 0, 0; m^ A ,m ix ,m ix ) . 



(74) 
(75) 



(76) 
(77) 



(78) 
(79) 
(80) 
(81) 



(82) 
(83) 
(84) 
(85) 



Note that the dipole form factors (74-77) are proportional to a fermion mass. The chirality 
flip takes place in the external fermion lines, for the terms proportional to LL and RR mixings 
and in the internal fermion lines (charginos or neutralinos), for the terms proportional to the 
LR mixing. 



The branching ratio £j — > reads 
BR(£j - 



a) = BR(£j -> tjvjv,) x w 

Lip 
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X 



AX ™% 



L[C]* L\C] T-i/ \ . x.a Liur n[u\ j-i i \ 

c iax CjaxF^xax) + c IAX c JA ^F 2 {x AX ) 



m xiL[cy 



R[C] 



L[N]* L[N] F , o 
C /AX c JAX r Z\ X A 



, m X A L[NYR[N] F(0 x 
^ ™ -TAX C JAX r 4:\ x AX) 

m ej 



+ R 



(86) 



References 

[1] M. L. Brooks et al. [MEGA Collaboration], Phys. Rev. Lett. 83, 1521 (1999) [hep-ex/9905013]. 

[2] K. W. Edwards et al. [CLEO Collaboration], Phys. Rev. D 55, 3919 (1997). 

[3] S. Ahmed et al. [CLEO Collaboration], Phys. Rev. D 61, 071101 (2000) [hep-ex/9910060]. 

[4] U. Bellgardt et al. [SINDRUM Collaboration], Nucl. Phys. B 299, 1 (1988). 

[5] D. W. Bliss et al. [CLEO Collaboration], Phys. Rev. D 57, 5903 (1998) [hep-ex/9712010]. 

[6] R. Akers et al. [OPAL Collaboration], Z. Phys. C 67, 555 (1995). 

[7] P. Abreu et al. [DELPHI Collaboration], Z. Phys. C 73, 243 (1997). 

[8] T. Yanagida, Prog. Theor. Phys. 64, 1103 (1980), also in Proceedings of the Workshop on the 
Unified Theory and the Baryon Number in the Universe, Tsukuba, Japan, February 13-14, 1979, 
eds. A. Sawada and A. Sugamoto, p. 95, KEK Report 79-18; M. Gell-Mann, P. Ramond and 
R. Slansky, in Proceedings of the Workshop on Supergravity, Stony Brook, New York, September 
27-28, 1979, eds. P. van Nieuwenhuizen and D. Friedman (North-Holland, Amsterdam, 1979), p. 
315-321; R. N. Mohapatra and G. Senjanovic, Phys. Rev. Lett. 44, 912 (1980). 

[9] J. I. Illana and T. Riemann, Phys. Rev. D 63, 053004 (2001) [hep-ph/0010193]; J. I. Illana, 
M. Jack and T. Riemann, hep-ph/0001273. 

[10] D. Atwood, S. Bar-Shalom, G. Eilam, and A. Soni, hep-ph/0203200. 

[11] J. A. Aguilar-Saavedra et al. [ECFA/DESY LC Physics Working Group Collaboration], "TESLA 
Technical Design Report Part III: Physics at an e + e~ Linear Collider", hep-ph/0106315. 

[12] G. Wilson, talks at DESY-ECFA LC Workshops in Frascati, Nov 1998 and Oxford, March 1999. 

[13] J. Hisano, hep-ph/0204100. 

[14] M. J. Levine, Phys. Rev. D 36, 1329 (1987). 

[15] M. Frank and H. Hamidian, Phys. Rev. D 54, 6790 (1996) [hep-ph/9603222]. 

[16] J. Hisano, T. Moroi, K. Tobe and M. Yamaguchi, Phys. Rev. D 53, 2442 (1996) [hep-ph/9510309]; 
J. Hisano and D. Nomura, Phys. Rev. D 59, 116005 (1999) [hep-ph/9810479]; S. w. Back, T. Goto, 
Y. Okada and K. i. Okumura, Phys. Rev. D 64, 095001 (2001) [hep-ph/0104146]. 

[17] D. F. Carvalho, M. E. Gomez and J. C. Romao, Phys. Rev. D 65, 093013 (2002) [hep-ph/0202054]. 



20 



[18] D. Delepine and F. Vissani, Phys. Lett. B 522, 95 (2001) [hep-ph/0106287]. 
[19] M. Frank, Phys. Rev. D 65, 033011 (2002). 

[20] I. Hinchliffe and F. E. Paige, Phys. Rev. D 63, 115006 (2001) [hep-ph/0010086]. 

[21] N. Arkani-Hamed, H. C. Cheng, J. L. Feng and L. J. Hall, Phys. Rev. Lett. 77, 1937 (1996) 
[hep-ph/9603431]; J. Hisano, M. M. Nojiri, Y. Shimizu and M. Tanaka, Phys. Rev. D 60, 055008 
(1999) [hep-ph/9808410]; J. Kalinowski, in Proc. of the APS/DPF/DPB Summer Study on the 
Future of Particle Physics (Snowmass 2001) ed. R. Davidson and C. Quigg, hep-ph/0202043. 

[22] A. Ilakovac and A. Pilaftsis, Nucl. Phys. B 437, 491 (1995) [hep-ph/9403398]. 

[23] J. Bernabeu and A. Pilaftsis, Phys. Lett. B 351, 235 (1995) [hep-ph/9502296]; A. Pilaftsis, Phys. 
Rev. D 52, 459 (1995) [hep-ph/9502330]. 

[24] P. Langacker and M. Plumacher, Phys. Rev. D 62, 013006 (2000) [hep-ph/0001204]. 

[25] E. O. Iltan and I. Turan, Phys. Rev. D 65, 013001 (2002) [hep-ph/0106068]. 

[26] C. x. Yue, H. Li, Y. m. Zhang and Y. Jia, Phys. Lett. B 536, 67 (2002) [hep-ph/0204153]. 

[27] F. Borzumati and A. Masiero, Phys. Rev. Lett. 57, 961 (1986); F. Gabbiani and A. Masiero, 
Nucl. Phys. B 322, 235 (1989); F. Gabbiani, E. Gabrielli, A. Masiero and L. Silvestrini, Nucl. 
Phys. B 477, 321 (1996) [hep-ph/9604387]. 

[28] M. Apollonio et al. [CHOOZ Collaboration], Phys. Lett. B 466, 415 (1999) [hep-ex/9907037]. 

[29] D. E. Groom et al. [Particle Data Group Collaboration], Eur. Phys. J. C 15, 1 (2000). 

[30] J. L. Feng, K. T. Matchev and Y. Shadmi, Nucl. Phys. B 613, 366 (2001) [hep-ph/0107182]. 

[31] D. F. Carvalho, J. R. Ellis, M. E. Gomez and S. Lola, Phys. Lett. B 515, 323 (2001) [hep- 
ph/0103256]; Z. Chacko and G. D. Kribs, Phys. Rev. D 64, 075015 (2001) [hep-ph/0104317]; 
F. Deppisch, H. Pas, A. Redelbach, R. Riickl and Y. Shimizu, hep-ph/0206122. 

[32] T. Moroi, Phys. Rev. D 53, 6565 (1996) [Erratum-ibid. D 56, 4424 (1996)] [hep-ph/95 12396]. 

[33] M. Knecht and A. Nyffeler, Phys. Rev. D 65, 073034 (2002) [hep-ph/0111058]; M. Knecht, 
A. Nyffeler, M. Perrottet and E. De Rafael, Phys. Rev. Lett. 88, 071802 (2002) [hep-ph/0111059]; 
M. Hayakawa and T. Kinoshita, hep-ph/0112102; I. Blokland, A. Czarnecki and K. Melnikov, 
Phys. Rev. Lett. 88, 071803 (2002) [hep-ph/0112117]; J. Bijnens, E. Pallante and J. Prades, 
Nucl. Phys. B 626, 410 (2002) [hep-ph/0112255]; J. Prades, hep-ph/0108192v3. 

[34] H. N. Brown et al. [Muon g-2 Collaboration], Phys. Rev. Lett. 86, 2227 (2001) [hep-ex/0102017]. 

[35] G. 't Hooft and M. Veltman, Nucl. Phys. B 44, 189 (1972). 

[36] G. Passarino and M. Veltman, Nucl. Phys. B 160, 151 (1979). 

[37] T. Hahn and M. Perez- Victoria, Comput. Phys. Commun. 118, 153 (1999) [hep-ph/9807565]. 

[38] G. J. van Oldenborgh and J. A. Vermaseren, Z. Phys. C 46, 425 (1990). 



21 



[39] T. Appelquist and J. Carazzone, Phys. Rev. D 11, 2856 (1975); A. Dobado, M. J. Herrero and 
S. Penaranda, Eur. Phys. J. C 12, 673 (2000) [hep-ph/990321lj. 

[40] W. Hollik, J. I. Plana, S. Rigolin, C. Schappacher and D. StSckinger, Nucl. Phys. B 551, 3 (1999) 
[Erratum-ibid. B 557, 407 (1999)] [hep-ph/9812298]. 

[41] G. W. Bennett [Muon g - 2 Collaboration], hep-ex/0208001. 

[42] K. Inami, T. Hokuue and T. Ohshima [Belle Collaboration], hep-ex/0210036. 



22 



